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LEGEND  FOR  FIGURES  OUTSIDE  THE  TEXT 


Fig.  3 


F(B)  versus  6.  (There  is  a  symmetry  with  respect  to  0  =  45°. )  Same  maximum 
value  *  1.59  at  6  ■  57'  for  1  MeV,  0  *  19*  for  9  MeV  and  6  «  11'  for  25  MeV. 


Same  angle  0  *  45°  for  mininum  value  *  4.52  x  10“^  for  1  MeV,  minimum 
value  •  2.46  x  10”^  for  9  MeV  and  minimum  value  *  7.51  x  10“  ^  for  25  MeV. 


Fig.  4 


(a) 


(b) 


F(0)  versus  0.  Same  maximum  value  *  1.59  at  0  *  57'  for  1  MeV, 

0  *  19'  for  9  MeV  and  0*11*  for  25  MeV.  Same  converging  value  * 

1  at  0  -  0°. 

F(0)  versus  0.  Similar  to  (a)  at  complementary  angles  (0  ♦  90°  -0). 


Fig.  6 


Ej  (impact  energy)  versus  E (beam  energy);  Ej  in  eV,  Eb  in  MeV;  with 
space  charge. 
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Fig-,,.? 

Ej  (impact  energy)  versus  Eb  (beam  energy);  Ej  in  eV,  Eb  in  MeV;  no  space 
charge. 
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Fig.  8 

Cross-section  for  ionization  of  atomic  hydrogen  taken  from  Kieffer  and  Dunn, 
Rev.  Mod.  Phys.  38,  1  (1966) 


Fig.  9 

No  space  charge;  beam  radius  (in  cm)  vs.  distance  (in  km)  for 
0.5  attenuation  for  various  propagation  angles. 
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Fig.  10 

No  space  charge;  beam  radius  (in  cm)  vs.  distance  (in  km)  for 
0.1  attenuation  for  various  propagation  angles. 
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Fig.  11 

No  space  charge;  beam  energy  (in  MeV)  vs.  distance  (in  km)  for 
0.5  attenuation  for  various  propagation  angles. 
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No  space  charge;  beam  energy  (in  MeV)  vs.  distance  (in  km)  for 
0.1  attenuation  for  various  propagation  angles. 
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Fig.  13 

Including  space  charge:  beam  radius  (in  cm)  vs.  distance  (in  km)  for 
0.5  attenuation  for  various  propagation  angles. 


(a)  IMeV;  1,000  mA/cm2 

(b)  IMeV ;  10,000  mA/cm2 

(c)  lOMeV ;  1,000  mA/cra2 

(d)  lOMeV ;  10,000  mA/cm2 

(e)  lOOMeV;  1,000  mA/cra2 

(f)  lOOMeV;  10,000  raA/cm2 

From  top  to  bottom:  5°(85°),  10#(80°),  45° 


Fig.  14 

Including  space  charge:  beam  radius  (in  cm)  vs.  distance  (in  km)  for 
0.1  attenuation  for  various  propagation  angles. 


(a)  IMeV;  1,000  mA/cm2 

(b)  IMeV;  1 0,000  mA/cra2 

(c)  lOMeV ;  1,000  mA/cra2 

(d)  lOMeV ;  10,000  mA/cm2 

(e)  lOOMeV;  1,000  mA/cra2 

(f)  lOOMeV;  10,000  mA/cra2 

From  top  to  bottom:  5°(85°),  10°(80°),  45° 


f.lfr.Jl 

Including  space  charge:  beam  energy  (in  MeV)  vs.  distance  (in  km)  for 
0.5  attenuation  for  various  propagation  angles. 

(a)  55  cm;  1,000  mA/cra2 

(b)  55  cm;  10,000  mA/cra2 

(c)  100  cm;  1,000  mA/cm2 

(d)  100  cm;  10,000  mA/cra2 

From  top  to  bottom:  5C(85°),  10°(80°),  45° 

(NOTE:  No  effect  is  discernable  for  beam  radius  of  10  cm) 


Fig.  16 

Including  space  charge:  beam  energy  (in  Mev)  vs.  distance  (in  km)  for 
0.1  attenuation  for  various  propagation  angles. 


(a)  55  cm;  1,000  mA/cra2 

(b)  55  cm;  10,000  mA/cm2 

(c)  100  cm;  1,000  mA/cm2 

(d)  100  cm;  10,000  mA/cra2 

From  top  to  bottom:  5#(83°),  10°(808),  45° 
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INTRODUCTION 


When  a  neutral  particle  beam  passes  through  the  upper  atmosphere,  it 
undergoes  a  number  of  physical  processes  which  lead  to  the  degradation  of  the 
beam  intensity  and  of  the  beam  energy.  Probably  the  most  important  amongst 
these  processes  is  stripping  which  leads  to  the  loss  of  (one  or  more) 
electrons  of  the  beam  particle.  Stripping  can  be  induced  by  collisions  with 
neutral  particle  ions  or  electrons,  without  substantial  differences  in  the 
stripping  cross  section  per  electron.  While  it  is  usually  assumed  that  the 
stripping  particles  originate  from  the  atmospheric  gas  or  plasma  through  which 
the  beam  passes,  this  is,  as  we  pointed  out  and  discussed  in  an  earlier 
Scientific  Report,  [l]  (to  be  referred  to  as  I)  not  necessarily  so.  The 
electrons,  stripped  off  the  beam  particles,  themselves  can  become  agents  of 
subsequent  stripping  processes.  Instrumental  in  this  scenario  is  the  presence 
of  the  earth's  geomagnetic  field  which  helps  to  accumulate  the  beam  generated 
electrons  in  the  beam’6  path.  It  was  actually  Alfven  [2]  who  realized  that 
the  passage  of  a  neutral  beam  through  magnetic  field  brings  about  an  anomalous 
source  of  Ionization. 

In  I,  certain  basic  aspects  of  this  beam  induced  stripping  (BIS)  were 
analyzed.  We  demonstrated  through  a  rather  simplified,  but  qualitatively 
faithful  model,  that  at  high  enough  beam  current  BIS  may  become  the  dominant 
mechanism  of  beam  degradation. 

Two  major  conditions  have  to  be  satisfied  in  order  for  the  BIS  process  to 
become  effective:  (i)  the  stripped  electrons  must  remain  within  the  beam  for 
a  long  enough  time  so  that  they  can  undergo  further  ionizalng  collisions,  and 
(11)  the  stripped  electrons  must  acquire  a  high  enough  (but  not  too  high) 
relative  velocity  with  respect  to  Che  beam  particles,  so  that  they  can  cause 
stripping  collisions. 
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The  presence  of  the  geomagnetic  field  is  critical  in  meeting  both  of 
these  conditions.  In  addition  to  the  direct  effect  of  the  magnetic  field,  the 
separation  of  the  positive  ions  and  negative  electrons  causes  a  polarization 
electric  field  which,  in  combination  with  the  magnetic  field  generates 
complicated  drift  trajectories. 

The  effect  of  the  magnetic  field  can  be  especially  clearly  seen  in  the 
case  of  the  beam  propagating  along  the  magnetic  field.  In  a  stripping 
collision,  the  electrons  are  ejected  at  a  large  angle  with  respect  to  the  beam 
velocity  and  thus,  without  the  magnetic  field,  would  promptly  leave  the  beam. 
The  magnetic  field,  however,  forces  the  electrons  back  into  the  beam.  For  a 
general  beam  orientation  the  situation  is  more  complicated;  nevertheless,  tne 
general  tendency  of  the  electron  becoming  trapped  by  the  magnetic  (or  by  the 
combined  magnetic  and  electric)  fields,  still  prevails.  In  the  other  extreme 
situation,  when  the  beam  propagates  perpendicular  to  the  magnetic  field,  it  is 
the  polarization  field  that  forces  the  electrons  to  travel,  for  a  substantial 
period  of  time,  with  the  beam.  This  scenario,  although  somewhat  similar  to 
the  well-known  problem  of  a  plasmoid  or  plasma  beam  travelling  across  the 
magnetic  field  [ 3]  is  different  in  that  the  electrons  are  continuously 
created  along  the  beam,  and  thus  the  critical  transient  phase  of  the  plasma 
beam  problem,  where,  upon  penetration,  the  beam  has  to  create  the  field  that 
propels  it  further,  plays  a  much  less  important  role. 

Having  demonstrated  in  I  that  the  BIS  process,  under  certain 
conditions,  can  be  an  important  factor  in  beam  degradation,  we  may  now  list 
the  principal  problems  that  one  has  to  address  in  order  to  obtain  a  detailed 
description  and  assessment  of  the  BIS  process. 


(a)  Calculation  of  the  “geometric  factor",  i.e.  the  fraction  of 
time  spent  by  the  electron  within  the  beam. 

(b)  Calculation  of  che  relative  velocity  of  the  electrons  and 
determination  of  the  relative  energy-dependent  stripping 
cross-section. 

(c)  Calculation  of  the  polarization  electric  field. 

This  Report  addresses  itself  to  item  (b)  and  to  one  aspect  of  item  (c). 

In  section  II  electron  trajectories  in  a  crossed  electric  and  magnetic  field 
configuration  are  calculated  without  the  usual  assumption  of  the  electric 
field  being  uniform  -  a  model  valid  only  in  the  case  of  small  gyro-radii.  It 
is  shown  that  with  strong  enough  spatial  variation  of  the  electric  field  the 
concept  of  ExB  drift  doesn't  apply  at  all.  In  Sections  III  and  IV  we  examine 
in  some  detail  the  generation  of  the  relative  velocity  by  the  magnetic  field 
and  find  that  it  depends  significantly  on  the  beam  energy,  on  the  beam 
diameter  and  on  the  propagation  angle  of  the  beam.  Then  we  calculate  the  beam 
decay  distance  for  10%  and  50%  beam  degradations,  as  a  function  of  beam  energy 
and  propagation  angle.  In  all  these  calculations  we  use  two  scenarios: 

(i)  no  polarization  electric  field  and  (ii)  a  model  electric  field 
representing,  in  the  absence  of  more  elaborate  self-consistent  calculations, 
the  polarization  electric  field;  for  the  geometric  factor,  again  in  the 
absence  of  a  detailed  calculation,  we  assume  the  median  value  y  *  0.5.  In 
Section  V,  we  give  some  qualitative  discussion  of  the  mechanism  determining 
the  actual  value  of  y. 
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II.  TRAJECTORY  EQUATIONS 

The  mathematical  problem  of  Lorentz  force  equation  coupled  with  Poisson 
equation  cannot  be  solved  analytically.  A  self-consistent  field  approach 
involve  making  appropriate  physical  approximations  and  numerical  calculations. 
However,  under  certain  non-trivial  conditions,  the  Lorentz  force  equation  with 
uniform  magnetic  field  and  external  electric  field  can  be  solved  exactly.  The 
resulting  electron  and  proton  trajectory  equations  are  useful  In  providing 
estimates  for  the  beam-induced  stripping  (or  self-stripping.) 

We  will  study  in  detail  two  cases:  (1)  beam  propagation,  uniform  magnetic 
field  and  linear  electric  field  are  mutually  orthogonal;  (2)  bean  propagation, 
uniform  magnetic  field  and  constant  electric  field  are  arbitrarily  oriented. 
These  two  cases  are  exactly  solvable  and  their  solutions  provide  useful 
insights  for  future  work. 

(1)  Mutually  Othogonal  v^,  B  and  E(y) 

*  >  ♦ 

We  designate  vb,  B  and  E(y)  respectively  as  the  beam  propagation  velocity, 

geomagnetic  field  and  linearly  varying  electric  field.  Specifically,  we  have 
♦  A  +  A  ♦ 

vjj  *  v^i ,  B  *■  Bk  and  E(y)  ■  (E  -  ey)j,  where  E  and  e  are  constants  to  be 
specified  later. 

The  Lorentz  force  equation  for  electron  becomes 
••  • 
x  “  -wy 


where 


•  eE  .  et 

mx - +  —  y 

m  m 


eB 


(1) 


<d  ■  —  ,  and  m  is  electron  mass. 

m 


Thi9  is  an  exactly  solvable  problem  and  the  resulting  electron  trajectory 
equations  are: 


(t)  e_]l/2  sinflt  +  (vb-u)REe)t 


(2) 


y(t)  -  Rge  (l-cosnt)  +  y0  (3) 

where  the  modified  frequency  0  is: 

Q  -  u»  (1  -  ~g) »  (4) 

and  the  modified  radius  Rge  is 


t  [yb  -  j  (E~ey0>) 

REe  *  l-l-  *  (5) 

u)B 

We  have  assumed  the  boundary  conditions:  x(0)  ■  0 
and  y(0)  »  y0,  |y0]  <  R0,  where  R0  is  the  beam  radius. 

Eqs.  (4)  and  (5)  indicate  the  effect  of  linearly  varying  electric  field 
with  a  non-zero  e.  Reduction  to  the  simple  case  of  uniform  field  can  be  seen 

c  B  1  E 

by  setting  e  *  0,  and  we  have  Q-hj  ■  —  and  Rge  ♦  —  [vj,  -  — ] . 

♦ 

The  Lorentz  force  for  proton  with  an  external  linear  field  E(y)  *  (E+ay) 
becomes 


y 


.  eE 
u)px  ^  jj 


y 


(6) 


where  Wp  ■  M  »  proton  mass,  and  the  constants  E  and  a  are  to  be 


specified  later.  The  resulting  proton  trajectory  equations  are: 


xp(t) 


RE 

f  “Vll/Z  *  8inV  *  (vb-»pREa>t 

1  "  u.pB-1 


(7) 


yp(t)  -  Rea  (cosflpt  -  1)  +  y0 


(8) 


where  the  modified  frequency  and  radius  are  respectively 
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A  simple  model  of  uniformly  distributed  space  charges  (electrons  and  protons) 
resulting  from  beam-atmosphere  ionizations  yeild  the  following  expressions  for 

E(y), 


E(y)  -  E  -  ey 


en 


£o^bR0 


2Re  +  Ro 
Re  +  "Rq 
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Re  +  Ro 


y]  . 


(id 


R0  <  y  <  2Re  +  R0 


E(y )  -  E  -  ey 


(12) 


en 
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e04LR0  LRe  +  R0 
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en 
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(13) 


-  (2Rp  -  R0)  <  y  <  R0 

where  Rp  and  Re  are  respectively  proton  and  electron  gyroradius,  L  is  about 
6Rp,  and  n  is  the  total  number  of  ionizations  for  beam  traversal  distance  L. 

The  results  obtained  so  far  are  based  on  simple  assumptions,  and  are  not 
adequate  for  physical  applications.  However,  they  constitute  a  basis  for 
further  considerations. 


£ 
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(2)  Arbitrarily  Oriented  vb,  B  and  Constant  E 

*■  * 

As  shown  in  Fig.  1,  we  have:  vb  ■  vbi, 


Figure  1 


written  as: 

x  -  -Yx  “  wzy 

y  =  -Yy  -  u,x2  +  U )2X 

2  -  -Yz  +  wxy 

i  v  -  v  *  e^v  y  c  z 

where  Yx  -  Yy  -  Yz  m  , 

eBY  eBz 

“x  *  V*  and  ^  "  ~V  * 


B  *  Bxi  +  Bzk 


E  -  Exi  +  Evj  +  Ezk 


and  tan  6 


The  Lorentz  force  equation 
equation  for  electron  can  be 


This  system  of  coupled  equation  can  be  solved  exactly  and  the  resulting 
electron  velocity  equations  are: 


x(t)  »  ~<*»ZA  sin  (<ot  +  $) - 7  (<oxYx  + 

U\*- 


<o_Yv  »Z? 

+  — 2*  +  (1  -  — )vb 


z(t)  *  toxA  sin  (tot  +  4>)  “  ^J7  (wxYx  +  uzYz)t 


where 


+  ^x 

2  2  »> 

0)^  0)^ 


y(t)  ■  toA  cos  (tot  +  $)  +  “2  (“xYz  “  “z^x^* 
A  -  1  [(Yy  -  <ozvb)2  +  2|]  1/2 1 

IJA  ' 


tan*  -  jj-  (wzvb  -  Yy)» 
u2  ■  +  Uy2  an<j  o  .  u»xyz  -  uizyx  • 

The  electron  trajectory  equations  can  be  obtained  by  integrating  the  above 
equations  with  the  boundary  conditions:  x(0)  •  0,  y(0)  -  y0,  z(0)  -  z0,  with  yQ 
^  R0  and  Zq  ^  Rq • 

The  Lorentz  force  equation  for  proton  can  be  expressed  as: 

x  -  -yx  “  “>zy 

y  ■  -*Yy  -  d»xz  +  cjxx  (18) 

Z  -  -yz  +  <*>xy 


where 


'Lx  «y  ecz 

Tx  *  — •  Yy  "  M  *  Yz  *  — ‘ 


ux  ”  m  «  ant*  uz 


Tne  proton  velocity  equations  are  given  by: 


9  ^  X 

x(t)  *  +uzA  sin  (ait  +  *)  +  — j  (<*)xYx  +  (*>zYz)t 
4^  4  (1  ^f2  )»/ 


(t)  *  -WjjA  sin  (wt  +  $)  +  -y  (u>xyx  + 

(il  __V  tA.JiS  _ 


U)xYy  , 

- +  — -  vb 


y ( t )  -  uiA  coz  (<ut  +  *)  +  (u)xYz  •  wzyx), 

uiz 


where 


*■  iT[(Yy-„Jvb)2  4£i]  1/2  . 


tan*  -  j  (ujzvb  -  Yy), 

o)2  m  ^^2  +  ^^2  an(j  d  .  ujxyz  -  o)zYx  • 

Again  the  proton  trajectory  equations  can  be  easily  obtained  by  integrating 


the  above  equations  with  the  same  boundary  conditions  as  those  for  the  electron. 


•  • 

It  is  worth  noting  that  the  accelerating  terms  in  x(t)  and  z(t)  for  both 
electron  and  proton  contain  the  factor  (“>XTX  +  <*>zyz),  which  can  be  considered  to 
have  resulted  from  the  presence  of  Ex  and  Ez.  The  sinusoidal  and  constant  terms 
represent  steady  state  motion  with  bounded  energy.  Therefore,  we  require  that 
(<*>xyx  +  wzYz)  vanishes  with  non-vanishing  Ex  and  Ez.  This  can  be  accomplished  by 
making  the  physically  reasonable  assumption  that  for  self-consistent  and 
self-sustaining  motion  there  exists  only  transverse  "polarizing"  electric  field, 
i  .  e .  , 


The  proton  velocity  equations  are  rewritten  as: 
xp(t)  *  -tan©  ( 8in  (“pt> 

-  sin0  -sin8.v^  cos  (uipt) 

+  sin  6.  +  cos20.(vb) 

zp(t)  «  !*■  sin  (oipt) 

+  cos0  -  sin0*vb^  cos(a)pt) 

-  cos0  (^j  +  co80*sin0  (vb) 

yP(t)  ■  ll~co8(v)J 

+  ~sin0.vb^  sin(uipt) 

where  <*>p  -  (<*>x2  +  <*>z2)i/2  -  J  (Bx2  +  Bz2)1/2  -  jjp 


(25) 


(26) 


(27) 


It  is  worth  noting  that  both  the  electron  and  proton  velocities  along 
direction  are  equal  to  cos0.vb.  More  specifically  we  have 
•  • 

(xe>n  +  (ze)«  "  co80.vb  -  V|( 

and 

(Xp)||  +  (z  p )  |  -  co80.vb  -  Vj| 

where  v>  is  component  of  beam  veloc'  y  vb  along  the  B  direction. 
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III.  BEAM-INDUCED  STRIPPING 


We  designated  the  relative  electron-beam  velocity  as  veb  -  ve  -  vb,  whose 
components  vex,  vey  and  vez  are  given  by  Eqs.  (22)  to  (24)  as 


E* 

'ex  “  *e  "  vb  "  tan®  g”  8in(u»et) 


-  sin8* 


|vbsln0  -  ME  -  cos(u>pt)J 


(28) 


'ey 


‘  Bliss’  [*  ‘  C0S("»t:i] 


+  £vbsin0  -  sin(tuet) 


(29) 


vez  *  -  b“  8ln(u,et> 


+  cos9  |vbsin0  -  -  cos(u>et)j 


«►  ♦  * 


(30) 


For  the  relative  proton-beam  velocity  vpb  ■  vp  -  vb,  the  components  are 
given  by  Eqs.  (25)  to  (27)  as 

Ez 

vpx  •  vp  -  vb  -  -  tan0*  sinUpt) 


Evl  r 

-  sin0*  vbsln0  -  1  -  cos(a)pt) 


(31) 


'py 


■  b7os9*  [*  ' 
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(32) 


~^vj,8in0  “  sin(uipt) 


'pz 


(r-) 


+  cos©*  (vfcSinG  -  jl  -  cos(<Dpt)|  (33) 

The  beam  (neutral)  particles  all  travel  with  velocity  vjj  within  the  beam 

2 

confine  which  is  of  cylindrical  form  with  cross-sectional  area  xR  .  The 


o 

charged  particles  (electrons  and  protons)  move  with  velocities  whose 
components  are  given  by  Eqs.  (22)  to  (27).  However,  the  charged  particles 
spend  only  a  short  time  within  the  beam  confine.  We  designate  the  time 
interval  as  Ty  and  consider  the  kinematics  of  the  charged  particles  within  the 
beam  confine.  As  shown  in  Fig.  2(a), 


a  charged  particle  starting  from 
the  beam  center  (i.e.,  origin  of 
x-y-z  system)  can  be  regarded  as 
staying  within  the  bean  confine 
if  its  moving  center's  path  is 


<  Du  - 


*  sin©' 


As  shown  in  Fig.  2(b),  vy  is  the 

♦ 

component  of  beam  velocity  v^,  in 
the  B  direction,  i.e., 

Vy  *  VfcCOsB.  We  regard  Ty  as 
the  time  interval  during  which 
the  charged  particle  is  within 
the  beam  confine.  Thus  we  have 


'i 


V|  v^alnB'cosO  vbsin20 


We  now  assume  Chat  all  Che  charged  particles  starting  from  the 
2 

cross-sectional  area  «Rq  (at  x  ■  0)  spend  approximately  the  same  time  interval 
T|  within  the  beam  confine.  For  this  time  interval  Tf ,  the  time-averaged 
relative  electron-beam  and  proton-beam  velocity  components  as  given  by  Eq. 


(34)  to  (39)  become 


<vex>  ■  tanO1 


h.  -L  Ti . 
B  *e  L 


cos(*e) 


-  sin0*  vbsin0  -  — 


ir.-i. 

JL 


sin(i^e) 


<v_v>  *  _ —  • 

e”  BcosO 


f-t 


sin( 


.5* 


+  vbsln0  -  - 


]k  [- 


cos(i(ie) 


<vez>  “  -  f  *  cos(*e>J 


+  COS0* 


[vbSin0  -  jl  -  ^  sin(^e>j 


<VpX>  ■  -  tan0 


•  r  ^  -  c°9<*p)] 


.ft 


-  sinG •  vb8in0  -  - 


[l '  ^  “"‘v] 


2eBR0 

where  \ Jv  -  u)eT||  -  - rr;  and 

e  e  '  mvbsin20 


(38) 


(39) 


2eBRo 

♦p  wpT|i  Mvjj8in20 


We  expect  electron-beam  collision  to  occur  for  <|vejj|>  *0,  i.e.  ,  non-zero 
time-averaged  magnitude  of  the  relative  electron-beam  velocity  vector 


veb  "  ve“vb*  from  Eqs.  (36)  to  (40),  we  have 


< I vebl >  "  [<vex>2  +  <vey>2  +  <vex>2]  1/2  “  <veb> 


where  the  right-hand  side  can  be  simplified  in  the  following  manner  to  give  an 

♦ 

upper  limiting  value  for  < | ve b I > * 

Firstly,  we  let  Ey  ■  v^BsinO  »  v^Bx,  i.e.  ,  the  electric  field  Ey  is  such 
that  the  charged  particles  move  with  the  beam  with  velocity  v^.  In  a 
self-consistent  space-charged  polarization  field,  this  would  represent  the 
upper  limiting  value  for  Ey.  Secondly,  we  let  Ez  ■  v^Bcos©,  which  represents 
the  optimum  Ez  for  charged  particles  having  drift  velocity  vj,  in  y-direction, 
which,  however,  is  balanced  by  Ex  ■  -  v^BsinO,  on  account  of  tranverse  field 
condition  EXBX  +  EZBZ  ■  0. 


Similarly,  proton-beam  collision  occurs  for 

♦ 

^lvpb!^  *  0»  From  Eqs.  (34)  to  (39),  we  can  obtain 
<veb>2  *  <lvebl>2  “  <vex>2  +  <vey>2  +  <vez>2 

and  <vpb>2  -  <|vpb|>2  -  <vpx>2  +  <vpy>2  +  <vpz>2 

We  will  now  consider  two  limiting  cases:  namely,  "transient"  case 
with  zero  electric  field,  and  "steady  state"  case  with  optimum  electric 
field. 

We  let  Ex  »  Ey  ■  E2  *  0,  and  Eqs.  (34)  to  (36)  become, 

<vex>o  "  vb  8ln2e  l(^T>  ain  «»e>  ~l] 

<vey>o  *  vb  sin  8  (^-  )  [ 1  -  cos  (*e)J 
<vez>o  "  vb  8ln  e«cos0  [(^-)  sin  (*e)  -l] 

and  from  which  we  finally  get 

-  sin20  {2  (~)2  [  1-cos  (*e>]  +  1-2  (^-)  sin  (*e)>  (40) 


where  the  subscript  "o"  denotes  quantities  with  zero  electric  field. 


Examination  of  Eq.  (40)  shows  that  the  factor 

{2(^-)2  [l-cos(*e)J  +  1_2(~— )  sin  (*e)}  -  F(0 ) 

'•'e  “e 

is  symmetric  with  respect  to  8  ■  45°  (also  the  angle  for  minima),  increases 
from  8  <  45°  (and  0  >  45°)  with  maxima  at  8  <  1'  (and  0  >  89°)  and  then 
oscillates  with  decreasing  amplitudes  and  increasing  frequency  and  finally 
converges  to  unite  at  0  ■  0®  (and  90°).  These  features  are  graphically  given 
by  Figure  3,  Figure  4(a)  and  Figure  4(b). 

The  8in20  term  in  Eq.  (40)  destroys  the  symmetry  of  F(6)  with  respect 
to  F(8)  with  respect  to  8  ■  43®.  The  quantity  <veb>20/vb2  approaches  zero  as 
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cos(weT)  -  cos(oe)  -  1  - 


Ro _ 

P.e  8in0 


Incorporation  of  the  above  results  leads  to  two  different  expressions 
for  <veb>2/v£  as  follows: 


-<-7b>°-  -  sin20  {2(f-)2  ll-cos(*e)]  +  l-2(-p-)  sin(*e)j 
v  fe  Vq 


where  = 


2R, 


e  Re  sin20 


with  6  <  90° 


(40) 


and 


^veb^o  *  sin20  •  (— )2  [l-cos(oe)]2 
- 5 -  0- 


(41) 


1  Kn 

where  oe  *=  cos-1  [1 - - - ] ,  with  0  =  90°. 

K-  sint) 


Similarly,  for  the  proton  case,  we  have 


£^b^o_  *  sin20  {2(y- )2  [  l-cos('J'p)  ]  +  l-2(~)  sin(^p)}  (42) 


2Rr 


Where  *P  *  Rp  Vin20 


1883 


-,  with  0  <  90' 


and 


— 0  -  -  sin20  •  (-i— )2  [l-cos(oD)]2 


(43) 
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where  oD  -  cos-1  (1  -  — - ■) ,  with  0  «  90°. 

v  Rp  sin0 

These  two  pair  of  equations,  Eqs.  (40),  (41)  and  Eqs.  (42),  (43), 
represent  the  results  for  the  "transient’*  case  where  the  electric  field  is 
zero,  i.e.,  beam  propagation  time  duration  is  too  short  for  the  formation  of 
space-charge  field. 

We  now  consider  the  "steady  state”  case  where  the  space-charge 
(polarization)  field  has  been  fully  extablished.  To  approximate  the  proper 
space-charge  field  which  is  to  be  obtained  from  a  self-consistency 
calculation,  we  propose  the  following  field  specifications: 

Ey  ■  vbB  sin0 

Ez  »  vbB  cos20  •  sin6 

Ex  *  -vbB  cos9  •  sin20 

These  field  components  agree  with  the  reasonable  physical  requirements 
as  described  in  the  following: 


(1) 

h  "  EZ  "  Ex  - 

U  for  0-0° 

(2) 

Ey  -  vbB,  Eg  - 

Ex  -  0  for  0  -  90° 

(3) 

'  >  i 

E :  “  vbB  s in0 

♦  | 

/  l+Cos20»  wit*1  maximum  E 

-  vjjB  at  0  -  90°. 

In  conditions  (2)  and  (3),  we  have  assumed  optimum  polarization  field 

strength  so  that  Ey  (and  E  )  ■  vbB,  i.e.,  the  charged  particle  velocity  is 

■¥ 

equal  to  beam  velocity  vl.  Substituting  these  field  components  in  Eqs.  (34) 


to  (39),  we  obtain  the  corresponding  "steady  state"  results  as  follows: 


— 7— -  -  -^(sin20)2{2(-j—  )2  [  1  —cos ( «^e >  ]  +  1~2(^— )  sin^e)}  (44) 

v,  ^  *e  ’'e 

D 

2R 

Where  +*  “  "R7sin20  *  Wlth  6  <  90°» 

— — “4-(sln28)2  (— )2  [l-cos(ce)j2  (45) 

v.  4  °e 

D 

R 

where  oe  *  cos"*  [1 - -  °  "s"n ~) »  with  6  ~  90°, 

^  2 

— jjr^ -  •  4-  (sin28  )2  (2(-j— )2  ll-cos(^p)  ]  +  l-2(-—)  sin(yp)}  (46) 

vfa  H  *p 

2RD 

where  *P  "  ~lTp'~sTn20’  With  6  <  9°° 


and 


<y>2 

Vb 

where  op 


t  (sin28)2  (~—  )2  (l-cos(op)]- 
°P 

:  c°s"‘  11  ’  -TpTile  ’•  "lth  8  *  90°- 


(47) 


To  obtain  the  electron-beam  "impact"  energy  E^0  for  the  transient 
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field-free  case,  we  first  let  H(oi,8,f)  -  <veb>0  /v^  as  given  by  Eq.  (40)  and 
write 


eb>£ 


j  mv^  H(a , 0 , f ) 
*o 

10cm 


(48) 


where  a 


and  f 


Similarly,  the  electron-beam  impact  energy  Ej  for  the  steady-state 
space-charge  field  case  can  be  written  as 

Ex  *  y  m<veb>2  -  y  mvb2  G(a,0,f)  (49) 

where  G(a,0,f)  is  given  by  Eq.  (44)  with  the  same  as  given  in  Eq.  (48).  Ej 
can  again  be  expressed  as  a  function  of  Eb,  R0  and  0.  It  should  be  mentioned 
that  Ej  is  symmetric  with  respect  to  6  *  45°. 

Figures  6(a)  to  6(d)  show  Ej  versus  Eb  (with  9even  values  of  R0)  for 
q  _  2° ,  5°,  10°  and  45°  respectively.  Figures  7(a)  to  7(q)  show  Ejq  versus  Eb 
(with  eleven  values  of  R0)  for  6  =  5°,  10°,  15°,  . ..,  85°  respectively.  The 
important  results  from  Figs.  6(c)  and  6(d)  can  be  stated  as  follows: 

(1)  The  critical  beam  radius  R0  is  55cm,  with  an  impact  energy  of  16e 

and  an  ionization  cross-section  Oj  equal  to  0.1  (na2).  [Note  that 

ita2  ■  1.13673  (10~16  Cm2)].  These  values  are  constant  for  Eb 
o  J 


ranging  from  IMev  to  100  Mev. 

(2)  For  6  ranging  from  10°  to  80°  and  Eb  from  5  Mev  to  100  Mev,  R0  and 
the  corresponding  ionization  cross  section  Oj  are  tabulated  in  the 
f ol lowing 


Ro  in  cm 

60 

70 

80 

90 

100 

oT  in  u a2 

1  0 

j  0.2 

0.5 

0.6 

0.7 

0.8 

(3)  For  0  <  5°  (and  0  >  85°),  Ej  decreases  for  E^  <  40  Mev.  The 
decreases  are  more  pronounced  for  larger  R0,  smaller  0  and  Eb, 
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The  important  features  shown  by  Figs.  7(a)  to  7(r)  are  as  follows: 

(1)  For  the  same  values  of  beam  radius  R0,  the  impact  energy  Ejq 


increases  with  the  angle  6. 

(2)  For  the  Rc  and  6,  Ej0  is  almost  constant  for  beam  energy  E^  ranging 
from  10  MeV  to  100  MeV. 

(3)  For  Efc  ranging  from  1  MeV  to  10  MeV,  Ej0  decreases  with  E^,  and  the 
decreases  are  pronounced  from  larger  R0. 

(4)  For  6  *  5°,  the  critical  R0  i6  again  55cm  for  E^> 1 0  MeV,  and  for 
6  -  80°,  E10  is  18  eV  for  R0  -  10cm. 

Because  of  large  proton  to  electron  mass  ratio  M/m  =  1843,  we  have  noil 

2 

result  for  the  case  of  ionization  by  protons.  In  other  words,  <Vpb>o/VbZ  and 
<vpb>2/vb2  as  given  by  Eq.  (42)  and  Eq.  (46)  respectively  are  less  than  8xl0~® 
(even  for  the  largest  R0  ■  100cm),  with  resulting  impact  energies  less  than  the 
threshold  ionization  energies. 
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IV.  ESTIMATE  OF  THE  GEOMETRIC  FACTOR 


In  this  Report  our  principal  aim  has  been  to  elucidate  the  dependence  of 
the  degradation  of  the  electron  impact  energy  on  the  different  beam  parameters. 
However,  as  discussed  in  the  Introduction,  the  BIS  process  requires  that  two 
conditions  be  met:  the  electrons  have  to  acquire  sufficient  energy  for 
ionization  and  they  have  to  be  retained  for  a  sufficiently  long  time  within  the 
beam  to  be  able  to  interact  with  the  beam  particles.  To  be  more  precise,  the 
electrons,  in  general,  leave  and  re-enter  the  beam  region  many  times.  Thus,  in 
I  we  have  defined  a  "geometric  factor"  Y,  Y  <  1  which  describes 
phenomenologically  the  reduction  of  the  collision  probability  due  to  the 
reduced  time  the  electrons  spend  within  the  beam.  The  calculation  of  the 
geometric  factor  for  a  general  situation  will  be  the  subject  of  another  Report. 
Some  comments,  however,  are  in  order. 

The  mean  free  path  of  an  electron  with  respect  to  ionizing  collisions  with 
the  beam  particles  is  obviously  quite  long,  comparable  to  the  beam  degradation 
length.  For  beam  propagation  at  very  small  angle  the  magnetic  field  can  force 
the  electrons  to  re-enter  the  beam  many  times  (even  though  they  acquire  a 
perpendicular  velocity  component  in  the  stripping  process)  and  to  effectively 
travel  with  the  beam.  Details  of  this  process  were  discussed  in  I.  For  an 
arbitrary  propagation  angle,  it  is  the  combined  effect  of  polarization  electric 
field  and  of  magnetic  field  that  is  expected  to  cause  the  electron  to  re-enter 
and  travel  with  the  beam.  The  details  of  this  scenario  are  to  be  discussed  in 
a  separate  Report.  However,  for  the  purpose  of  the  calculations  in  the  present 
Report  we  have  ignored  the  re-entry  of  the  electrons  and  calculated  averages 
over  the  short  time  while  the  electron  in  which  the  beam  after  its  emergence. 

We  believe  that  because  of  the  quasiperiodicity  of  the  motion,  this  doesn't 
seriously  affect  the  results  concerning  the  velocity  averaging. 
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It  should  be  noted  that  the  conditions  relating  to  the  relative  velocity, 
on  the  one  hand,  and  to  the  geometric  factor,  on  the  other,  are  in  general, 
contradictory.  Extended  stay  within  the  beam  requires  traveling  with  the  beam 
particles  and  thus  low  average  relative  velocity.  The  difficulty  is  partly 
due,  however,  to  the  model  used  for  our  calculation,  where  the  relative 
velocity  is  identified  with  the  relative  average  velocity.  Thus,  in  extreme 
situations,  the  latter  may  become  zero,  while  the  former  does  not.  Obviously, 
a  more  rigorous  calculation  is  called  for  in  such  cases. 

The  most  difficult  problem  in  a  rigorous  handling  of  the  particle 
trajectories  would  be  the  self-consistent  determination  of  the  polarization 
electric  field.  Since  such  calculation  is  much  beyond  the  scope  of  the  present 
Report,  in  Section  III  we  have  used  a  simple  physical  model  for  the 
representation  of  the  electric  field.  This  model  includes  a  component  parallel 
to  the  beam  velocity  (always  perpendicular,  however,  to  the  magnetic  field). 

The  justification  for  the  existence  of  such  a  parallel  component,  which  is  not 
contemplated  in  the  usual  theories  of  cross-field  propagation  [3],  can  be  found 
through  a  simple  glance  at  the  geometry  generated  by  a  neutral  particle  beam 
moving  at  an  angle  to  the  magnetic  field  and  undergoing  stripping.  Over  a 
distance  small  compared  to  the  ion  gyro-radius,  the  beam  carries  a  net  positive 
charge  (since  the  ions  deviate  little  from  their  original  path),  while  the 
electrons  create  a  triangular  wedge,  fanning  out  along  the  beam.  It  is 
expected  that  such  a  geometry  will  indeed  give  rise  to  a  parallel  component 
of  the  electric  field. 
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V.  CONCLUSIONS 

The  following  conclusions  emerge  from  this  study: 

1.  As  already  discussed  in  I,  the  BIS  process  can  become,  at  high 
altitudes,  a  non-negl igible  contribution  to  beam  degradations.  While  in  I 
we  discussed  the  case  of  propagation  parallel  to  the  magnetic  field,  in 
the  present  Report  we  considered  a  general  propagation  direction  at  an 
arbitrary  angle  to  the  magnetic  field. 

2.  The  propagation  of  the  beam  at  an  angle  to  the  magnetic  field  causes 
the  electrons  generated  by  the  BIS  process  to  acquire  sufficient  relative 
velocity  to  induce  further  stripping.  The  stripping  cross  section  by 
electrons  being,  however,  a  sensitive  function  of  the  impact  energy  (see 
Fig.  8),  which  in  turn,  is  a  function  of  the  beam  energy,  angle  of 
propagation  and  beam  diameter,  the  beam  attenuation  coefficient  becomes  a 
function  of  all  these  parameters,  in  addition  to  the  beam  current  density 
which  determines  the  density  of  the  beam  particles.  These  dependences  are 
shown  in  Figs.  9  through  16.  There  are  two  major  kinematic  approximations 
used  in  these  calculations:  (a)  representing  the  relative  velocity  by  the 
average  relative  velocity  and  (b)  calculating  averages  over  the  first 
portion  of  the  particle's  trajectory  whthin  the  beam.  Apart  from  extreme 
situations,  both  of  these  approximations  are  reasonable  and  also  can  be 
easily  removed  in  future  work. 

3.  The  polarization  electric  field  plays  an  important  role  in  guiding  the 
electrons  along  the  beam.  The  determination  of  this  field  for  an 
aribtrary  propagation  angle  is  an  extremely  complicated  task:  no  attempt 
has  been  made  in  this  direction  in  the  present  Report.  We  have  U6ed, 
however,  a  reasonable  model  for  representing  the  polarization  field  and  we 
have  calculated  effects  both  without  (corresponding  to  a  "transient" 
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situation)  and  with  polarization  fields.  The  differences  can  be  assessed 
by  comparing  Figs.  3  and  4,  Figs.  6  and  7  and  Figs  9-12  and  13-16, 
respectively.  In  Section  II. 1  we  have  also  examined  the  effect  of  the 


inhomogeneity  of  the  electric  field  and  have  found  that  if  the  scale 

E  1 


length  of  the  inhomogeneity,  say  Y  satisfies  Y  <  g-  —  the  usual  drift, 
trajectories  change  into  open  trajectories  and  lead  to  an  unbounded  motion 
which  cannot  be  described  in  terms  of  the  customary  drift  trajectories. 
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